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ABSTRACT 

We investigate the possibility of accounting for the currently inferred 
primordial abundances of D, ^He, ^He, and ^Li by big bang nucleosynthesis in 
the presence of baryon density inhomogeneities plus the effects of late-decaying 
massive particles (X), and we explore the allowed range of baryonic fraction 
of the closure density fib in such context. We find that, depending on the 
parameters of this composite model (characteristic size and density contrast 
of the inhomogeneities; mass-density, lifetime, and effective baryon number in 
the decay of the X particles), values as high as flbh^Q — 0.25 — 0.35 could be 
compatible with the primordial abundances of the light nuclides. We include 
diffusion of neutrons and protons at all stages, and we consider the contribution 
of the X particles to the energy density, the entropy production by their decay, 
the possibility that the X-products could photodissociate the light nuclei 
produced during the previous stages of nucleosynthesis, and also the possibility 
that the decay products of the X-particles would include a substantial fraction 
of hadrons. Specific predictions for the primordial abundance of Be are made. 

Subject headings: cosmology: theory — dark matter — early universe — nuclear 
reactions, nucleosynthesis, abundances 
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1. Introduction 

Standard homogeneous big bang nucleosynthesis could have produced the 
observationally inferred primordial abundances of D, ^He, ^He, and ^Li, provided 
that the baryon fraction of the cosmic closure density flf, would lie in the range: 

0M<n,hl<0.08 (1) 

where h^Q is the Hubble constant in units of 50 km s~^ Mpc~^ (Walker et al. 1991; Copi, 
Schramm, & Turner 1995). For the long-time most favored cosmological model, a flat 
Universe with Qm = 1 and Qa = (those being, respectively, the fractional contributions 
of matter and vacuum energy densities to the closure density), the upper bound to Qj, 
would mean that most matter in the Universe should be in nonbaryonic form. Given the 
far-reaching implications of the dominance of nonbaryonic dark matter, possible alternatives 
to homogeneous big bang nucleosynthesis have been explored, especially during the last 15 
years or so. 

The suggestion that the quark-hadron phase transition might be first-order and 
generate baryon inhomogeneities (Witten 1984) led to the calculation of the possible effects 
on primordial nucleosynthesis (Applegate & Hogan 1985; Applegate, Hogan, & Scherrer 
1987; Malaney & Fowler 1988). The goal was to see whether inhomogeneous big bang 
nucleosynthesis with Qh = 1 might account for the primordial light-element abundances. 
Besides a first-order quark-hadron phase transition, other mechanisms might also generate 
baryon inhomogeneities. Much of the work in this line is reviewed by Malaney & Mathews 
(1993). However, the recent studies, treating accurately the coupling between baryon 
diffusion and nucleosynthesis, show that the upper limit on flh set by the light-element 
abundances does not significantly differ from that obtained for homogeneous big bang 
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nucleosynthesis (Mathews, Schramm, & Meyer 1993; Thomas et al. 1994). This last 
conclusion, though, has very recently been challenged by Orito et al. (1997), who explore 
the dependence of primordial nucleosynthesis on the geometry of baryon inhomogeneities 
and find that cylindrical geometry might allow to satisfy the observational constraints with 
baryon fractions as high as ^bh^o ^0.2. 

A different approach has been to explore the possible modifications of the yields from 
homogeneous big bang nucleosynthesis by the effects of the decay of unstable massive 
particles (M few GeV) , produced at earlier stages in the evolution of the Universe and 
with half-lives longer than the standard nucleosynthesis epoch {tx >, 10^ s) (Audouze, 
Lindley, & Silk 1985; Dommguez-Tenreiro 1987; Yepes & Dommguez-Tenreiro 1988; 
Dimopoulos et al. 1988). Gravitinos produced during reheating at the end of inflation are 
a possible example of such particles. In Dimopoulos et al. (1988), the emphasis is put on 
the resulting hadron cascade. The main problem encountered in this model is the predicted 
overproduction of ^Li: ^Lif^Li » 1, whereas observations show ^Lif^ Li < 0.1. 

Although they have only been considered separately, baryon inhomogeneities and the 
presence of unstable massive particles decaying when the Universe has already cooled down 
below Tg ~ 0.4 are by no means mutually exclusive. Here we explore their combined effects 
on the primordial abundances of the light elements. The parameter space now has, of 
course, a dimension which is the sum of those for the two separate cases: characteristic size 
and density contrast of the inhomogeneities, mass-density, lifetime, and mode of decay of 
the massive particles. We find that there are regions in such extended parameter space 
where values of Jl;, as high as ^bh^Q — 0.35 would still be compatible with the primordial 
abundances of the light nuclides inferred from observations. Such values of fib are of the 
same order as the low values for Qm currently derived from a variety of sources, including 
high-redshift supernova searches (Perlmutter et al. 1998; Garnavich et al. 1998). Our 
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results thus suggest again the possibihty that all the matter in the Universe could be 
baryonic. 

On the other hand, recent determinations of the D abundance in high-redshift QSO 
absorbers, when confronted with the currently inferred primordial ^He abundance, might 
be in conflict with the predictions of standard, homogeneous big bang nucleosynthesis 
for Ni, — 3 (Steigman 1998): the "low" high-redshift D abundances (which appear more 
reliable) would indicate too high a value of Qf, to be compatible with that corresponding to 
the ^Hc abundance. Since it is hard to tell whether this conflict points to new physics or just 
to systematic errors in the derivation of abundances, Steigman, Hata, & Felten (1998) have 
discarded the constraint on from standard big bang nucleosynthesis and turned to other 
observational constraints to determine the key cosmological parameters. The results from 
our composite model, by showing how minor deviations from the standard hypotheses can 
produce agreement with the primordial abundances inferred from observations, support that 
attitude. Besides, as we will see, the combined effects of inhomogeneities plus late-decaying 
particles might solve the conflict between D and ''He abundances. 



2. Model, Results, and Discussion 

In the present model, nucleosynthesis first takes place in a Universe with baryon 
inhomogeneities. Later, when the temperatures are low enough for the chemical abundances 
to be frozen, unstable massive particles start to decay producing both hadronic and 
electromagnetic showers which alter the abundances of the light elements resulting from the 
previous stage. 

For the inhomogeneities, we consider a simple model consisting of two types of zones, 
one with high density and the other with low density, characterized by their respective 
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volume fractions /„ and 1 — fy, and by the density contrast R between the two zones 
(see Thomas et al. 1993). The inhomogeneities, produced at some earher stage, lead to 
differential diffusion of protons and neutrons when the temperature, in the expanding 
Universe, drops below T ~ 1 MeV and protons and neutrons are no longer in equilibrium 
(Applegate & Hogan 1985). Then, due to the much longer mean free path of the neutrons, 
the initial baryon inhomogeneities transform into variations of the local neutron/proton ratio 
(Applegate, Hogan, & Scherrer 1987, 1988; Rauscher et al. 1994). When nucleosynthesis 
starts at Tg ~ 1, neutrons are already uniformly distributed whereas the protons retain 
the spatial distribution they had at weak decoupling. Nucleosynthesis thus takes place in 
two different types of zones: the proton-rich and the neutron-rich ones. Since almost all 
neutrons end up into ^He, with the rapid growth of the abundance of this nuclide protons 
become exhausted in the neutron-rich zones and the same occurs with neutrons in the 
proton-rich zones. A density contrast thus appears again and neutrons start to diffuse back 
from the neutron-rich zones into the proton-rich ones. We have coupled neutron diffusion 
with the nuclear reaction network as in Rauscher et al. (1994). 

On the other hand, we have the X-particles, with masses and half-lives t^- The 
particles are massive {nix > 10 GeV) and we consider only half-lives in the interval 
10'^ s < Tx < s. The former means that they are nonrelativistic well before the start of 
nucleosynthesis. The lower limit to the half-life implies that the thermonuclear reaction 
rates have dropped to zero before the X-particles start to disintegrate, while the upper 
limit ensures that their decays leave no imprint on the cosmic background radiation. Prior 
to decaying, the X-particlcs just give a contribution to the matter-energy density. We 
define r as the number ratio of the X-particles to photons at some fiducial temperature Tq. 
Here we follow the history of cosmic expansion starting at Tq — 10^^ K. The product: 



is one of the parameters of the model. When a particle of mass higher than a few GeV 
decays during the keV era, a fraction of its energy goes into electromagnetic decay products 
while another fraction goes into hadrons. Both kinds of products give rise to showers 
as they thermalize with the background plasma. The baryonic decay products interact 
with the ambient protons and a-particles and initiate chains of nuclear reactions that 
lead to production of the hght nuclides D, ^He, ^Li, ^Li, and to destruction of some ^He 
(Dimopoulos et al. 1988). The electromagnetic decays consist of the injection of high-energy 
photons, electrons, and positrons into the background plasma. The high-energy photons 
dissociate the light nuclei. The average number of baryons ub produced in the decay of 
the X-particles has been calculated by Schwitters (1983). The numbers of nuclei i and 
neutrons produced by the disintegration of each X-particle can be calculated by modeling 
the hadronic cascades (Dimopoulos et al. 1988, 1989). 

In order to calculate the final abundance for each light nuchde i, including the effects of 
inhomogeneities plus X-particle decay, we integrate, starting at Tq = 10^^ K, the following 
set of differential equations: 



— ^ = ±SinKYn H -exp{-t/T^) 

at Tx 

where the + sign corresponds to the proton-rich zone, the - sign to the neutron-rich 
one, K is the diffusion rate of neutrons, rjo is the initial n^/n-y ratio, is the effective 
baryonic branching ratio of the X-decays, t^{E) is the photon spectrum produced by 
disintegration of the X-particles, a^i{E) is the photodissociation cross-section, and ac{E) 
is the Coulomb scattering cross-section. The integrals extend from the photodissociation 



ac(E) 



dE 



+ 



dt 



(3) 



ftf.d 
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threshold energy Qi to E^nax — f^t/i'^^T) (which thus increases with time). Subscripts 
n and p refer to neutrons and protons, respectively, and the subscript std indicates the 
standard contribution of thermonuclear reactions to dYi/dt. In order to deal with diffusion 
in the way indicated in (3), the nuclear abundances Yi = Ni/Ni, are taken relative to the 
total number of baryons in the given volume before neutron diffusion sets in, as in Rauscher 
et al. (1994). The neutron diffusion rate is given by: 

n = ^^Tl'\l + Q.7lQT,f^ (4) 
[d/a) 

where d/a (in cm MeV) is the comoving length scale of the inhomogeneities (Applegate, 
Hogan, & Scherrer 1988). The effective baryon branching ratio r*^ in (3) is another 
parameter of the model. It takes into account the dependence of the number of baryons 
produced by disintegration of the X-particles on their mass mx together with the 
dependence of the yields on the kinetic energies of the primary shower baryons: 

r% = (y) tbF (5) 

where is the true baryonic branching ratio, and F incorporates the dependence of the 
yields on the kinetic energy of the primary shower baryons (our "standard" have been 
calculated for nix — 1 TeV, vb — 5, and = 5 GeV). The photons produced by the 

decay of the X-particles have a spectrum: 



Our model, therefore, has five extra free parameters in addition to r]o'- the volume 
fraction fy and the density contrast R characterizing the inhomogeneities (but we also 
vary the comoving length scale d/a governing neutron diffusion), the mass m^, (in fact the 
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product rm-r) and the half-life of the X-particles, and the effective baryon ratio in 
the X-particle decays. We thus solve the equations (5) for different combinations of those 
parameters, searching for the regions in the parameter space where the predicted primordial 
abundances of the light nuclides are compatible with those inferred from observations. The 
nuclear reaction network consists of 161 thermonuclear reactions, plus 20 reactions involved 
in the hadron cascades, plus 15 photodisintcgration reactions. In following the expansion of 
the Universe we take into account the contribution of the X-particles to the total energy 
density and the entropy change due to their disintegration (assuming that the X-decay 
products thermalize on a time scale much shorter than the expansion time scale). The 
cross-section sources and a more detailed description of the model will be given elsewhere 
(see also Lopez-Suarez 1997). 

The set of primordial abundances to be fitted is the following: 




(7a) 






(7c) 






(7e) 
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The limits in (7a-e) are adopted, respectively, from McCuUough (1992), Geiss & 
Reeves (1972), Pagel & Kazlaukas (1992), Krauss & Kernan (1994), and Smith, Lambert, 
& Nissen (1996). In (7e) we have adopted a conservative upper limit, taking account of 
the possibility that any primordial ^Li might have been partially destroyed (^Li remaining 
almost intact) in metal-poor halo stars. 

We have varied the volume fraction within the interval 0.01 < < 0.28, and the 
density contrast within 50 < i? < 5000. For the comoving length scale d/a, we have tried the 

cases without diffusion and with 10^'^ cm MeV < d/a < 10'^'^ cm MeV. For the parameters 
related to the X-particles, 10"^ GeV < rm^ < 10^ GeV, 1.5 x lO'^^ < < 1.5 x 10"^ 
and 10^ s<T^< 10^ s. 

As an example of the results, in Figure 1 we show the dependence of the final D 
abundance on Tx, for three different values of /„ and fixed values oi R, d/a, r'^, rm^, and 77/ 
(the final value of the baryon to photon ratio). We see that for 5 x 10^ s <Tx <,Q x 10^ s the 
resulting abundances fall within the range allowed by observations. In Figure 2 we show the 
dependence on of the final abundances of D, ^He, ^He, ^Li, ^Li, and ^Be, for fixed value 
of /„ and the same values of the other parameters as in Figure 1. The two vertical dashed 
lines mark the interval of values of Tx compatible with all the observationally inferred 
abundances of the light nuclides. 

In summary, wc obtain results which are compatible with the observations for little 
diffusion {d/a — 10^'^ cm MeV), small abundances of the X-particles {rmx ~ 10~^ GeV), 
and modest numbers and energies of the shower baryons (1.5 x 10~^^ < r% < 1.5 x 10~^^). 
The density contrast between the two model zones must be 500 < R < 5000, and 
the volume fraction 0.144 < /y < 0.192. The half-lives of the X-particles are 
6.19 X 10^ s < Ta; < 7.43 X 10^ s, as illustrated in Figures 1 and 2. In that region of 
the parameter space 18 < 7710 < 22 (7710 being here, as ususal, rj in units of lO"^'^). That 
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corresponds to a baryon fraction: 



0.25 < Qbhlo < 0.35 



(8) 



in strong contrast with (1). A testable prediction of the present model is the production of 

an appreciable amount of ^Be: 



The production of ^Be is characteristic of inhomogeneous big bang models (Malaney 
& Fowler 1989; Jedamzik et al. 1994; Orito et al. 1997). The current observational upper 
hmit to the Be abundance (Duncan et al. 1997; Garcia Lopez et al. 1998) is of the same 
order as the vahics found here. There is, however, the problem that B/Be ~ 10 almost 
down to [Fe/H] — —3, and since our model predicts B/Be < 1, agreement would require a 
drop in the ratio taking place below some still smaller metallicity. 



We have shown, by means of a simple model, that the combined effects on big bang 
nucleosynthesis of baryon inhomogeneities plus the decay of unstable, relatively long-lived 
massive particles, giving rise to both electromagnetic and hadron cascades, might be to 
allow agreement with the primordial light-element abundances inferred from observations 
for values of Qfc much higher that those allowed by standard, homogeneous nucleosynthesis. 
The upper limit might be as high as ^bh^Q — 0.35. The values obtained here are of the same 
order as the low values now being derived from a variety of sources and, therefore, they 
pose in new terms the question of whether all matter in the Universe could be baryonic. A 




(9) 



3. 



Conclusions 
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testable prediction of the model is the production of a ^Be abundance that is of the order 
of current observational upper limits. 

On the other hand, in the parameter region of our model where there is agreement 
between predicted and observationally inferred primordial light-element abundances, given 
values of flh (or, equivalently, tjiq) always predict "low" D abundances (in the sense of the 
high-redshift abundances referred to in the Introduction), thus potentially eliminating the 
conflict with the ^He abundance pointed out by Steigman (1998). 

The model presented here deals with inhomogeneities in a very simplified way. A 
futher step will be to examine the effects of the geometry of the density fluctuations on 
the outcome. Orito et al. (1997) have already shown that cylindrical shell geometry alone 
(without the extra effects of late-decaying particles) might allow fib ^ 0.2 (but for density 
contrasts R ~ 10^, much higher than those considered here). Another extension of the 
model will be to consider particles with shorter half-lives, decaying at the time when 
thermonuclear reactions are still taking place. 
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Fig. 1. — Primordial D abundance as a function of the t^, the half-hfe of the X-particles, 
for three different values of the volume fraction f^: 0.074 (short-dashed line), 0.119 (long- 
dashed line), and 0.144 (dot-dashed line), and fixed values of the other parameters (see text 
for the meaning of the different symbols). 
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Fig. 2. — Primordial abundances of the light nuclides D, ^He, ^He, ^Li, ^Li, and ^Be, as 
a function of the half-life of the X-particles, t^, for fixed values of the other parameters: 
ri^^^ j = 18, = 0.144, R = 5000, d/a = 10^•^ rm^ = 10'^ and r*^ = 1.5 x 10"^^ (see text for 
the meaning of the different symbols). The two vertical dashed lines mark the boundaries 
of the interval where the predicted abundances are compatible with those observationally 
inferred. 



